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Oh. Abstract 

. Group theoretical methods are used to study some properties of the Riccati equation, which 

is the only differential equation admitting a nonlinear superposition principle. The Wei-Norman 
method is applied for obtaining the associated differential equation in the group SL(2, R). The 
superposition principle for first order differential equation systems and Lie-Scheffers theorem are 
also analysed from this group theoretical perspective. Finally, the theory is applied in the solution 
of second order differential equations like time independent Schrodinger equation. 

o 

oo 1 Introduction 

On ' 

Nonlinear phenomena are having everyday more and more importance in almost all branches 
of science, and in particular in Physics. One of the most important nonlinear differential 
equation is the Riccati equation (see e.g. [TJ and references therein). This differential 
equation has recently been shown to be related with the factorization method (see e.g. 0, [|, 
f|, H). The recent interest of this equation is steadly increasing since Witten's introduction 
of supersymmetric Quantum Mechanics |6|, [7j]. 

Two important features of this Riccati type differential equation are: 
i) Even if we cannot find the general solution by means of a finite number of quadratures 
over elementary functions of the coefficients di(t) defining the differential equation, 

dx(t) 2 



Oh 



X 



dt 



a Q (t) + ai(t) x(t) + a 2 (t) x 2 (t) } (1.1) 



once a particular solution X\ (t) is known, the change of variable x = X\ + z leads to a new 
differential equation of the Bernouilli (for n = 2) type: 

dz 

— = (2a 2 Xi + a ± )z + a 2 z 2 . (1.2) 
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This is a particular instance of (|1.1|) for which a = 0. Notice that under the change of 
variable w = 1/x the Riccati equation ( |1.1| ) becomes 



dw{t) 



+ a {t) w 2 (t) + ai(t) w(t) + a 2 {t) = 0. 



(1.3) 



dt 



In particular, the Bernouilli equation ( |1.2|) can be reduced to a linear one just by intro- 
ducing the new variable u = 1/z. In this way, if we know a particular solution the general 
solution can be found through two quadratures. 

ii) When three particular solutions, Xi(t), x 2 (t), x 3 (t), of the differential equation (|1.1| ) 
are known, the general solution can be written with no quadrature at all: 



where k is an arbitrary constant determining each particular solution. 

In this sense we can say that there exists a nonlinear superposition principle for the Riccati 
equation, because the general solution can be expressed as a function x = Q(x\,X2,X3,k) of 
three particular solutions and one arbitrary constant k. 

Our aim is to explain these facts from a group theoretic viewpoint and present some new 
ideas both about the Riccati equation itself and on the nonlinear superposition principle, for 
which the Riccati equation is the simplest case. 

The organization of this paper is as follows. In Section 2 we review the problem of 
reducing a second order linear differential equation to a nonlinear first order Riccati equation, 
what means that the original linear superposition principle for the second order equation 
should be replaced by a nonlinear superposition principle. We also remark that this fact 
is due to a relation of the Riccati equation with the SL(2, M.) group to be explicited later. 
In Section 3 we explain a method developed by Wei and Norman || for determining the 
solution of a differential equation in a Lie group and we apply the method for the study of 
the Riccati equation, finding in this way the explicit form of the superposition principle as 
a consequence of some group theoretical computations. The superposition principle for first 
order differential equation systems and Lie-Scheffers theorem are studied in Section 4. It is 
shown that for an important class of such systems the problem of finding the general solution 
is reduced to the simpler problem of finding one particular solution of another system on a 
Lie group G, and moreover, even without solving directly this new system the solution of the 
original system can be easily found as soon as we know a fundamental set of solutions of it. 
Once again the simplest case is Riccati equation and the superposition principle can be found 
by determining the first integral of a system. Finally in Section 5 we give as an example the 
application of the Wei-Norman method in the solution of second order differential equations 
taking as a propotype the Schrodinger equation for the harmonic oscillator. 

2 The nonlinear superposition principle 

There is a well known method of relating a linear second order differential equation with a 
Riccati equation. Actually, given the linear second order differential equation 



Xy X> ]_ X> 2 Xl 




X — X2 X3 — X2 



d 2 u 

dp 



+ b(t) — + c(t) u = 0, 



(2.1) 
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the property of linearity means that the vector field 

x = «± 

Ou 

generates a one-parameter Lie group of point symmetries of the equation ( |2.1|) 

t(e) = t, u(e) = e e u. 

Changing coordinate u to a new one v = ip(u) such that the vector field X = ud/du 
becomes a translation generator (Straightening-out Theorem), i.e. X = d/dv, and therefore 
determined by the equation Xv = 1, leads to v = \ogu, i.e. u = e v , and then, 

du „ dv dv , 
— = e v — = u—. 2.2 
dt dt dt y ' 



When written in terms of this new coordinate the equation (|2.1p becomes 



d 2 v , , . dv f dv \ 2 

The unknown function v does not appear in the preceding equation and therefore a 
lowering of the order is obtained when introducing the change of variable x — and then 
we will get a Riccati equation 

dt 

as it was pointed out in |§. 

Therefore, the linear superposition principle for solutions of ( |2.1| ) translates in a nonlinear 
superposition principle for those of this Riccati equation, as it will be shown later. 

Notice that (|2~^ ) shows that 

1 du , n JS 

x=-—. 2.4 

u dt K ' 



The second order differential equation (|2 . 1| ) is equivalent to the set of fl2.3|) and 
We should also remark that two solutions ui,«2 of ( |2.1| ) project on the same solution of 
( p.3|) if and only if there exists a nonzero real number A G K such that ^2(0) = Ami(O) and 
^(0) = A«i(0). 

From the geometric viewpoint, Riccati equation can be interpreted as the one determining 
the integral curves of a time-dependent vector field 

d 

T = (a (t) + ai(t)x + a 2 {t)x 2 ) — . 

ox 

Let us remark that this vector field can be written as a linear combination with time 
dependent coefficients of the vector fields 

d d d 

L = —, L 1 = x— , L 2 = x 2 —, (2.5) 
ox ox ox 

which generate a three dimensional real Lie algebra with defining relations 

[Lq, Li] = L , [L ,L 2 ]=2L 1 [L x , L 2 ) — L 2 , (2.6) 
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and therefore isomorphic to sl(2, R). In fact, it is an easy matter to check that the (local) 
one-parameter transformation Lie groups of R generated by Lq, Li and L2 are 



x 

x 1— > x + e, inei, 37 1 — > 



\— xe 



i.e., they are fundamental vector fields corresponding to the action of <SX(2, R) on the real 
line R extended with a point at the infinity, R, given by 

•yx + \7 J 7 

and 

$(A, 00) = — , — ) = 00. 

7 7 

3 The Wei— Norman method 

Let us consider a differential equation system 

tfe*(t) 



dt 



X l (x,t), i = l,...,n, (3.1) 



which can be seen as the differential equation system whose solutions are the integral curves 
of the time-dependent vector field 

X = X i (x,t)^-. (3.2) 

The theorem for existence and uniqueness of solutions of the preceding differential equa- 
tion tells us that, for small enough t, there exists a map $ 4 applying the initial condition 
x(0) = (x l (0)) into the corresponding value x l (t). Correspondingly, functions / defined in a 
neigborhood of x(0) transform as 

[U(t)f](x) = f &\x)) , (3-3) 
and taking derivatives with respect to t we obtain 

d ( , . i \ / \ \ dx l d f 



dU{t) 
dt f 



« ((/ o •-)(,)) = ^(» r >(,)) 

*(/) (*r'w) = wiraw- (3.4) 



This relation is valid for any different iable function /, and therefore 

dU(t) 



dt 



U(t)X. 



We recall || that given such a differential equation for operators for X being a linear 
combination of vector fields in R n , 

m 
k=l 
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namely, 

^ = jr ak (t)U(t)L k , (3.5) 

where the L k span a finite dimensional real Lie algebra, it is possible to write the general 
solution in the form 

m 

U(t) = l[exp(g t (t)L t ), (3.6) 
i=i 

where the functions gi(t) are given by the solution of the system obtained from the relation 

Li, (3.7) 



^2ai(t)Li = J^dM 

i=l i=l 



II ^Pi-dji^adLj] 
j=i+i 



and the initial condition ^(0) =0, i = 1, ■ ■ ■ ,m. 

This method proposed by Wei-Norman (see also [ 10 1 ) can be used in the case of the Ric- 
cati equation and the generalization for other differential equation system involving several 
degrees of freedom is immediate. In fact, given the differential equation (|1.1| ) there will be an 
evolution operator U(t) which takes values in SL{2, R) and satisfies the differential equation 

= U(t)[a (t)L + a x (t)L x + a 2 (t)L 2 ], (3.8) 

together with the initial condition U(0) = I, and where the vector fields L k , for k = 0,1, 2, 
are fundamental vector fields associated to the left action of SX(2,K) on the extended real 
line R that are explicitly given in ( |2.5| ) and generate a Lie algebra isomorphic to sl(2, R). 

According to the Wei-Norman method ||, there will be functions u (t) , ui(t) , u 2 (t) such 
that 1*0(0) = Ui(0) = 1*2(0) = and 

U(t) = exp(wiLi) exp(u 2 L 2 ) exp(u L ). (3.9) 
Then, when using (|3.7f) , we obtain 



a (t)L + ai(t)Li + a 2 (t)L 2 = Ui(t) exp(— M (^)adL ) exp(— u 2 (t)g,dL 2 )Li 

+ u 2 (t) exp(-M (t)adL )L 2 + uo(t)L 

= ii\(Li — uqLq + u 2 L 2 — 2uqu 2 Li + UqU 2 L ) 

+ u 2 (L 2 - 2u L 1 + uILq) + u L , 

and therefore the system of differential equations for the functions Ui 

a o(t) = —UqUi + U 2 Q U 2 Ui + UqU 2 + ii 

ai(t) = iii — 2uqU 2 Ui — 2uqU 2 

a 2 (t) = u 2 u 1 + ii 2 , (3.10) 

that can be written in normal form, 

u Q (t) = a {t) + ai(t)ii (0 + a 2 {t)ul(t) 
ux{t) = ai(t) + 2a 2 (t)u (t) 

u 2 (t) = a 2 {t)-a l {t)u 2 {t)-2a 2 {t)u Q (t)u 2 {t). (3.11) 
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We remark that the first equation for uq is nothing but the original Riccati equation 
and therefore it seems that there is no advantage at all. However when looking in the other 
two equations, we see that provided that the appropriate solution for the first equation, 
the one determined by wo(0) = has been found, the solution for the second one is almost 
immediate and when introducing these values in the third equation this one becomes a first 
order differential equation and the solution reduces to a quadrature. In this sense we have 
reduced the problem of finding the general solution of (|1 . 1|) to the one of finding the particular 
solution such that u (0) = 0. This is quite similar to the property i) we mentioned in the 
Introduction. 

Once the solution of ( |3.11|) determined by u (0) = «i(0) = 1*2(0) = has been found, 



the general solution of the differential equation will be written 

x(t) = (U(t)x)(0) = [exp(MiLi)exp(M 2 ^2)exp(M Lo)x]| i= o, (3.12) 

and therefore 

<t)= eUlX ° ui +u , (3.13) 
1 — u 2 e Ul x 

where xq = x(0). 

Let us now fix three independent initial conditions. A possible set is the one given by 
£i(0) — ► 00, £2(0) = y £3(0) = 1, (actually any other three different initial conditions will 
be transformed to this one under an appropriate transformation of the group SL(2, M)). 

Having in mind the form ( |3.13| ) of the general solution, we see that the above mentioned 
initial conditions determine the particular solutions 

Xi(t) = ^— +u (t) 

x 2 (t) = u (t) 

gUl(t) 

xM = i-M»-"> + '" o(t) ' < 3 ' 14 ' 

and then the functions uo,Ui,u 2 are determined as 

u (t) = x 2 (t) (3.15) 
(x 3 (t) - x 2 (t))(x 2 (t) - Xi(t)) 



ui(t) = log 



(x 3 (t) - Xi(t)) 



(3.16) 



Mt) =~ ^tTWY (3 - 17) 

and therefore, when putting these values in fl3.13j ), we find that for the general solution 

= x xi(t)(x 3 (t) - x 2 {t)) + x 2 (t)(x 1 (t) - x 3 (t)) ^ lg , 

xo(x 3 (t) - x 2 (t)) + (x^t) - x 3 (t)) 



which is the well known superposition principle for Riccati equation ( |1 . 1|) which can also be 
written as 

(x - x 2 )(x 3 - x x ) 



(x - Xi)(x 3 - x 2 ) 



X . (3.19) 
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The factorization (|3.9| ) for the evolution operator is not the only possible one, but there 
are other five alternative factorizations. We next present the results for all possible reorder- 
ings: 

II) When we consider the factorization 

U(t) = exp(g L ) exp^Li) exp(g 2 L 2 ), 

the associated system turns out to be 

go = a e~ 91 
9i = ai- 2a g 2 
g 2 = a 2 - a x g 2 + a gj, 

with go(0) = f7i (0) = g 2 (0) = 0, and then the solution is 

e 91 {x + g ) 



(3.20) 



x(t) 



1 - g 2 e^ (x + g ) ' 
the three particular solutions X\,x 2 , x 3 being expressed as 

1 



(3.21) 



Xl(t) 

x 2 (t) 
x 3 (t) 



92 



9oe 



in 



1 - 9o92e 91 



o!Jl 



;i + go) 



with the inverse relation 



go 

9i 

92 



1- g 2 e^(l + g )' 

x 2 (xi - x 3 ) 
xi(x 3 - x 2 ) 

xj(x 3 - x 2 ) 



(3.22) 



log 



1 

Xi 



(xi - x 2 )(xi - x 3 ) 



(3.23) 



We should remark that the third equation in ( ^.20[ ) is a new Riccati equation and that 
once the solution for this new Riccati equation is found, we substitute its value in the second 
one and integrate without any difficulty, and when this value for g\ is put in the first equation 
we can integrate it. Therefore the new result we have found here is the following: if one 
knows a solution of a related Riccati equation, given by the third one in ( |3.20|) equation, 
satisfying g 2 (0) = 0, the general solution of ( |1 . 1| ) can be found in a straightforward way. 

Ill) The third possibility corresponds to the factorization 

U(t) = exp(h 2 L 2 ) exp(/iiLi) exp(h L ), (3.24) 
the associated system being 



ho = a + aiho + a 2 h\ 
hi = a\ + 2a 2 ho 
h 2 = a 2 e hl , 



(3.25) 
(3.26) 
(3.27) 
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and the general solution is then 



x(t) = - h h . 

1 - n 2 x 

The three particular solutions xi,x 2 ,x 3 will be written as 



(3.28) 



Xi(t) 

x 2 (t) 
x 3 (t) 



e hl 
h 2 

h 



1-h, 



+ h , 



(3.29) 



the inverse relations being 



h 


= x 2 




hi 


= log 




h 2 


x 2 - 


- x 3 




Xi - 


- x 3 



(x 3 - x 2 )(x 2 - Xi) 
(x 3 - Xi) 



Notice that in this approach the first equation in (|3.27|) is the original equation 
IV) The fourth reordering leads to the factorization 

U(t) = exp(JiLi) exp(/ £o) exp(/ 2 ^2), 



(3.30) 
O). 



(3.31) 



associated system 



fo — a o + a\fo — 2a /o/2 

/i = ai - 2a f 2 

f 2 = a 2 - a x f 2 + a f 2 , 



and general solution 



x(t) 



e fl x + fo 



l-f 2 {ehxo + fo) 
The three particular solutions are now 



(3.32) 
(3.33) 



xi(t) 
x 2 {t) 
x 3 (t) 



1 

J2 



fo 



1 - hh 



e h +fo 



l-/ 2 (e^+/ ) 



(3.34) 



and the inverse relation 



fo 
h 



x x x 2 
Xi - x 2 



log 



x\(x 3 - x 2 ) 
(xt - x 2 ){xi -x 3 ) 
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h = --• (3.35) 



We remark that now the third differential equation for f 2 in ( 3.32 ) is the same as in 
the case II) and it provides a new method of finding the general solution of ( |1 . 1| ) once the 
particular solution satisfying /2(C)) = of the associated Riccati equation is found. 

V) The fifth possibility is 

U(t) = exp(v L ) exp(v 2 L 2 ) exp^Li), (3.36) 

with associated system 

v = a Q e~ Vl 

v\ = a\— 2aov 2 e~ Vl 

v 2 = a 2 e Vl -a v 2 2 e' v \ (3.37) 



and solution 



1 - V 2 (x + Vq) 

the expressions for the three particular solutions xi, x 2 ,x% are 

xi(t) = 

v 2 

e Vl v 



x(t) = e " (X ; + "°\ ; (3.38) 



with the inverse relation 



v 



1 - V 2 V 

e^{l + v ) 
1 -v 2 {± +v ) 

(x 3 - xi)x 2 
(x 2 - x 3 )x 1 

x\{x 2 - x 3 ) 



vi = log 



(x 2 - x 1 )(xi - x 3 ) 



{x 2 -x 3 )xi 

v 2 = — -. (3.40) 

(x 2 - X 1 ){X 3 - Xi) 

VI) The last possibility is the factorization 

U(t) = exp(w 2 L 2 ) exp(w L ) exp(wiL 1 ), (3.41) 

and then the associated system is 

w = a e~ Wl — a 2 wle wi 
Wi = a\ + 2a 2 w e Wl 

w 2 = a 2 e w \ (3.42) 

the general solution 

x(t) = - °— + Wo e w \ (3.43) 

1 - w 2 x 
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and the expressions for xi, x 2 , x 3 , 



x 2 (t) 
x 3 (t) 



w 1 



W-2 



1 - w 2 



+ w e 



Wl 



+ w e 



U'l 



(3.44) 



with inverse relation 



w 

Wi 

w 2 



(x 3 - xi)x 2 
x 2 - xi)(x 3 - x 2 ) 
(x 2 -x t )(x a -x 2 ) 



log 



{x 3 - x 2 ) 
(x 3 - Xi 



(x 3 - Xi) 



(3.45) 



Finally, we remark that in these two last approaches there is no uncoupled differential 
equation of the Riccati type whose solution allows us to find the solution of the two other 
remaining equations in the system, and therefore the general solution of the original Riccati 
equation, anymore. However, if in the third equation in (|3.37|) we define <f>(t) = v\ we will 
get the Riccati equation 

' (3.46) 



-<p 2 + q(t)(f ) + p(t) 



with 



?(*) = — , 
a 



p(t) = di 



a i ■ a i 
2a a 2 a 2 + — . 

a 2 



(3.47) 



In a similar way, in the sixth case, taking derivatives in the second equation in ( p.42j ) 
and after some manipulations, the equation for cp — w\ becomes a Riccati equation 



where 



r(t) 



(p = s(t) + r(t)<p + - V 2 



Oi 2 . . d\ 

— , sit) = ai a 2 + 2a a 2 — 

a 2 a 2 



(3.48) 



(3.49) 



Let us summarize the results we have got. We have reduced the problem of finding the 
general solution of the Riccati equation to the one of determining a curve in SL(2, M) which 
is defined through its second class canonical coordinates, and this leads to a differential 
equation system. Once the curve in SL(2,M) is known we are able to find the general 
solution of the Riccati equation. However, even if we are not able to find the solution of the 
corresponding differential equation system for the second class coordinates in the group, we 
know the form (|3.12| ) of the general solution of the original Riccati equation. Even more, 
given a set of (three in the Riccati case) fundamental particular solutions we can determine 
the function giving us the superposition principle (|3.18|) . 

We have seen that the general solution of Riccati equation is given by 



x(t) 



x xi(x 3 - x 2 ) + x 2 {xx - x 3 ) 
x (x 3 -x 2 ) + (xi -x 3 ) 



(3.50) 
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where Xq is a constant depending on the initial conditions. We aim now to show how it 
is possible to reconstruct the original differential equation once the superposition formula is 
given. In the case of Riccati equation the superposition formula ( |3.50|) is equivalent to 

a — cx 

x ° = 1 1 

ax — o 

and it is easy to check that taking derivatives in the preceding relation we find 

(a — cx — cx)(dx — b) — (a — cx)(dx + dx — b) =0 

from where we find the following expression for x 

(cd — cd) 2 (—ad + ad + be — be) (ab — ab) 
(be — ad) (be — ad) (be — ad) ' 

that is a Riccati equation. 

On the other side, if we assume that there is a superposition formula for Riccati equation 
and we try to determine the function giving that formula, x = (j)(xi,X2, x 3 , k), we will have 

• d(j) . d(j) . d<f) . , j2 

x = (p = — — x± + - — x 2 + t; — %3 = clo + Oi0 + a 2 (p 

OXi OX 2 OX 3 

from where the following system of partial differential equations is found 

d(f) d(f) d<p 1 



dx\ 8x2 dxs 

dej) d(p d(p , 

Xlx hX 2 7 ^^3^— = 

OXi OX2 OX3 
50 90 dl _ 

Xl d Xl + 2 dx 2 + 3 dx 3 ~^- j 



Now a long computation leads to the following expression for 0: 



kxi(x 3 -x 2 ) +x 2 (x x -x 3 ) fo ko . 

x(t) = —-, — — — . (3.52) 

k(x 3 - x 2 ) + (xi - x 3 ) 



4 The superposition principle for first order differen- 
tial equation systems and Lie-Scheffers theorem 

We are now interested in studying the existence of a superposition principle for first order 
differential equation systems generalizing the one obtained for the Riccati equation. More 
explicitly, given a system 

==r = X>(x,t), i = l,...,n, (4.1) 
dt 

we ask whether there exists a set {x^ 1 ', . . . ,x^ mS) } of fundamental solutions and a function 
$ : R n ( m+1 ) _> R n such that the general solution of (O) can be expressed as 

x = $(x {1 \...,x {m) ;a u ...,a n ), 
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where Cbfi £1X6 constants related with the initial conditions. 

Before studying the general case, we remark that we know that, at least in the case where 
the system is the autonomous linear system 

dx . „ 
— = Ax, x e R ra , 
dt 

with A being a constant matrix, it defines a flow <p t which can be considered as a curve on 
the general linear group GL(WL,n) given by 4>t{x) = e At x. The flow satisfies 

d ^ A A 

and it can be determined from a fundamental set of solutions, i.e., when the n x n matrix 
X whose columns are the vectors defining the solutions {xi(t), . . . ,x n (t)}, 

X(t) = ( Xl (t),...,x n (t)), 

is an invertible matrix, then the equation 

X(t) = e tA X(0) 

shows that the evolution operator e tA is determined as e tA = X(t)X(0) _1 , i.e., the funda- 
mental system of solutions allows us to find the flow of our first order differential equation 
system. In other words, we have a possitive answer to our previous question with m = n 

and $ being the linear map . . . , x^; a± : . . . , a n ) = a\ x^ + h a n x^ n K 

From here it is clear that with any linear differential equation system on M. n we can 
associate an equation on GL(n, R) by setting 

9 = 9 A 

or 

9~ ld T = ^ 
dt 

The matrix A is an element of gl(n, R). 

Moreover, this way of finding the resolvent can also be used for time-dependent systems 

I = A{t)x - 

In this case if we have a fundamental set of solutions, denoted by X(t) = (xi(t), . . . ,x n (t)), 
from 

X(t) = R(t,0)X(0) 

we get R(t,0) =X(t)X(0)-\ 

We could equally well associate a time-dependent equation on GL(n, R) by setting 

with Ai being elements of the natural basis of gl(n, R). 
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The point we want to remark is that if the evolution preserves some structure, then 



A(t) = <f> t o — 

lies in certain subalgebra of gt(R, n). The flow now defines a one-parameter family of trans- 
formations, and conversely, given a one-parameter family of transformations it will determine 
a vector field in IR n , the corresponding fundamental vector field. 

The critical fact is that the general solution is determined by a linear operator from a 
set of fundamental solutions, and this is a linear superposition principle. It is very natural 
to ask what happens when the vector field is nonlinear. The answer is that, at least in some 
cases, to be explicited shortly, there is a kind of non-linear superposition principle, as it was 
proved by Lie This nonlinear superposition principle is simply a generalization of the 



previous construction to those cases where the action of the group is not linear and M. n is 
replaced for a manifold M. 

In the general case of the system ( |4.1| ), the Theorem for existence and uniqueness of 
solutions of such systems tells us that there will be, for each small enough t, a local dif- 
feomorphism of M n which establishes the correspondence among the initial values and the 
corresponding ones for the explicit value of the parameter t. In other words, the evolution is 
described by a curve gt in the group of diffeomorphisms of lR n . We have seen that when we 
consider a linear autonomous system this curve lies in the group GL(n,M) and is just the 
exponential of the matrix A giving the system, g(t) = exptA In the linear time-dependent 
case, we also have a curve g t in GL(n,M.) but it is not the exponential anymore: the only 
thing we can say is that 

,_i d<f> t 



A{t) 



1 dt 



takes values in the Lie algebra Ql(n, R). Actually, the solution is obtained by the Dyson, 
time-ordered, exponential. 

The point is that for other types of vector fields the curve described by gt belongs to 
other Lie subgroups of the group of diffeomorphisms and these are just the cases for which 
the idea of the nonlinear superposition principle is generalizable. 

The result established by Lie and Scheffers |TT| is that the general evolution defined by 



([Ll|) can be expressed in terms of m fundamental solutions if there are r vector fieds Yi, . . . , 
Y r , such that the vector field X, 

can be expressed as a linear combination 

X = a 1 (t)Y 1 + --- + a r {t)Y r (4.2) 
and furthermore the vector fields 

close a finite dimensional real (or complex) Lie algebra, with dimension r, i.e., there exist r 3 
real numbers c a p 7 such that 

[Y a ,Y p ]=c a ^Y 7 . (4.3) 

Moreover, in this case mn > r. When C a ( x ) = a% aj x ^ + Kx,-> with a l a j and b % a arbitrary 
constants, the system is linear. 
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Let us consider an effective action of a Lie group G of dimension r on n-dimensional 
differentiable manifold M, $ : G X M -> M, and $ 9 : M -> M and $ x : G -> M denote 
the maps = ^(sO = ^KflS^), for g G G, x G M. Choosing an initial point x(0), every 

curve g : I — > G in the group determine a curve in the manifold M by 

= $(<?(*), x(0)) = $ git) (x{0)) = % {0 )(g(t)), 

and taking derivatives with respect to t we see that the tangent vectors to the curves g(t) 
and x(t), respectively, are related by 

x(t) = $> x (o)* 9 (t)g(t). 

Let us remark that g(t) G T g (t)G and x{t) G T X ^M. 

We can express x(t) in terms of x{t): we recall that if X2 = & g (xi), then 

where R g denotes right translation in the Lie group G, because for any g' G G, 
$ X2 (g') = ®(g',x 2 ) = $(</, <&(<?, x x )) = ®{g'g,x x ) = o R g ){g'). 
Now, using the chain rule for computing the differentials we find that 

^X2* ±L g*i 

and then, 

We can use this relation for g = g(t) and x± = x(0) and then we find the following 
expression for x(t): 

X{t) = ®x(t)*e(Rg-i(t)*g(t)g(t))- 

Since R g -i^) is the right translation leading g(t) to the neutral element e G G and 
G T g ^G, then R g -i^ g u^g(t) G T e G and we know that T e G may be identified with the Lie 
algebra of G, q. Moreover, for linear Lie groups, i.e., subgroups of GL(n, R), right translation 
reduces to right multiplication by the corresponding matrix, and hence R g -^(t)*g{t)9if) is just 
the product g{t)g~ l {t). 

Let {e a } r a=1 be a basis of the corresponding Lie algebra g with defining relations 

and denote X a the corresponding fundamental vector fields defined by 

(X a f)(m) = — [/(exp(-te a )m)] |t=0 , 
for any differentiable function /. We recall that in this case 

[X a , Xp\ = C a /3 7 Xy. 

Now, if the time-dependent vector field defining the system ( }4.1|) is of the form 

X(t) = a a (t)X a 
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we associate with it the following differential equation on the Lie group itself 



g 1 ( t )~^ = a <x{t)e a - 

In this way the given system is replaced by a higher dimensional system of first order 
linear equations. Or in other words, we have replaced the original system of differential 
equations with a new system on the group G, but the important point is that it is enough 
to find a particular solution, the one starting from the neutral element, for obtaining the 
general solution of the system ( [Oj) : the solution starting from x is given by <&(g(t),xo). 

Moreover, even if we do not know the solution g(t) of the new system, it is possible to 
find it from the knowledge of a convenient set of particular solutions of ( |4.1|) such that ( |4.2| ) 
with the additional condition ( |4.3|) . More explicitly, given a curve X\{t) that is a particular 



solution of the given system, there are, in principle, different possible choices for the curve 
g{t) such that Xi(t) = $(g(t), Xi(0)), because the stability group of the point £i(0) may be 
nontrivial. If we choose a different particular solution, X2(t), then the ambiguity reduces to 
the group intersection of the isotopy groups of 2i(0) and £2(0). We can, if necessary iterate 
the procedure until we arrive to a set of m particular solutions xi, . . . ,x m allowing us the 
determination of the curve g{t). Of course as we have r unknown functions, the second class 
canonical coordinates, and we have mn conditions, it should be mn > r. 

More explicitly, a set of x\, . . . ,x m of solutions is said to be a fundamental system of 
solutions, if 

x 1 {t) = $(0(t),zi(O)) 

x m {t) = ®(g(t),x m (0)) (4.4) 

is a minimal set allowing us to solve for g(t) via the implicit function Theorem. If this can 
be done we get 

g(t) = F(xi(t), x m (t); £i(0), . . . , x m (0)), 
and then any other solution can be written as 

x(t) - $(F( Xl (t), . . . , x m (t); m(0), . . . , x m (0)),x(0)) = 0. 



Therefore the left hand side of this relation defines a constant of the motion. 

Starting with the action $ : G x M -*Mwe should find the minimal integer number 
m such that the isotopy group of the action of G on the product M m = M x ■ ■ • x M 
(m times), extended from $ by $ m (g, x x , . . . , x m ) = ($(g, Xi), . . . , x m )), reduces to the 
neutral element for a point such that any two coordinates are different. 

We recall that the fundamental vector field corresponding to an element of g generating 
a one-parameter Lie subgroup contained in the isotopy group of a point, vanishes in such a 
point, and conversely. Therefore, when expressed in terms of fundamental vector fields that 
means that the extensions to M m of the fundamental vector fields X a do not vanish in a a 
point whose coordinates are different. The general solution then is found by adding a new 
component and looking for constants of motion. 
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The procedure is next illustrated with an example for the simplest case n — 1. According 
to Lie's Theorem we should look for a finite dimensional real Lie algebra of differential 
operators 

X a = f a (x)^-. 

It can be shown that the only finite dimensional Lie algebra that can be found from 
vector fields in one real variable are sl(2,M) and its subalgebras. The uniquely defined (up 
to a change of variables) realization of sl(2, R) is given by 

2 & 

X Q = x—, X_ = — , X + = x —. 
ox ox ox 

These vector fields close the 5l(2,R) Lie algebra 

[X ,X.] = -X., [X ,X + ]=X + , [X_,X+] = 2X . 

Therefore, it suffices to consider the case in which the time-dependent vector field X 
defining the equation can be written as a linear combination X = a\X + a 2 X + + a X_ with 
°o — a>o(t), a>i = Q-i(t) and a 2 = a 2 (t) real functions. For simplicity we will consider the case 
when ai, a 2 and a 3 are real numbers and then we obtain the differential equation 

^ = X(x, t) — ao + a\x + a 2 x 2 , 

which is nothing but the well known Riccati equation. 
First we note that the determinant 



1 Xi x\ 



1 x 2 x\ 

4 



1 X3 x% 



(xi - x 2 )(x 2 - x 3 )(x 3 - Xi) 



is generically different from zero, while the system a + bx\ + cx\ = 0, a + bx 2 + cx\ = 
always has a solution. Consequently, m = 3 in this case. 

Now, for obtaining the general solution we should define the vector fields 

d d d d 
V o = X TT + X ^ ha; 3^— , 

OX OX\ OX 2 OX3 

d d d d 

and look for a solution of the system 

V f = V + f = V_f = 0. 

This system of partial differential equations is integrable, because the vector fields Vo, 
V + and V- close a Lie algebra, and therefore they define an integrable distribution. 

The last equation V_/ = tell us that the function / depends only on the differences 
ui — x — xi, u 2 — x\ — x 2 and u 3 = x 2 — x 3 , because the characteristic system is 

dx dx\ dx 2 dx 3 

T ~ T" ~ T" ~ "T' 
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and it has as first integrals the differences u-y — x — xi, u 2 — X\ — x 2 and u 3 = x 2 — x 3 . Now, 
if f(x,xi,x 2 ,x 3 ) = (p(ui, u 2 , u 3 ), the condition Vof = is written 

dip dip dip 

Ula ^ u ^a hM 3^— = 0, 

OUi OU2 OU3 

i.e., the function ip should be homogeneous of degree zero, and therefore it only can depend 
on the quotients v\ = u\ju 2 and v 2 = u 3 /u 2 , ip(u\,u 2) u 3 ) = (f)(vi,v 2 ). Finally, the condition 
V + f = can be written in these coordinates, after a long computation, as 

^ + 1)^-^ + 1)^ = 0. 
OV\ ov 2 

The corresponding characteristic system is 

dvi dv 2 



and taking into account that 



/ 



log 



we obtain that the constant of motion / should be a function of 

c - — V2 

Vi + 1 v 2 + 1 ' 

and therefore 

(x-x l )(x 2 -x 3 ) 

— c 

(x - x 2 )(xi - x 3 ) 

provides the non-linear evolution principle giving x(t) as a function of three independent 
solutions 

_ (xi - x 3 )x 2 c + xi(x 3 - x 2 ) 
(xi - x 3 )c + (x 3 - x 2 ) 

Let us remark that for n — 1 there is only one nonlinear differential equation family 
satisfying Lie-Scheffers theorem: the Riccati equation. Of course, proper subalgebras of 
s[(2, R) lead to the linear inhomogeneous equation when a 2 = or to a linear homogeneous 
equation when do = a 2 = 0. However, for n = 2 in addition to SL(3, M), 0(3, 1) and 0(2, 2), 
we can realize families of Lie algebras with arbitrary large Abelian ideals. 



5 Application in the solution of second order differen- 
tial equations 

Algebraic methods have very often been used in the search for eigenvalues of operators and 
the corresponding eigenvector. The particular case of the harmonic oscillator is the prototype 
and it is based on creation and annihilation operators, and therefore it is related with the 
Heisenberg group. The possiblity of relating linear second order differential equations with a 
Riccati equation, as indicated above, and the related group SL(2,M) has not been exploited 
till now, as far as we know. 



17 



In this section we will explore the use of Wei-Norman method based on the SL(2, R) 
group for studying the spectral problem of the second order differential operator determined 
for the Hamiltonian of the Harmonic oscillator 



H 



2M 



+ -X 2 , 
2 ' 



(5-1) 



where k is a constant. 

We will use the following notation: 

a - 

i -- 

A = 

and then the Hamiltonian can be written 

H 



k/M, 
ax, 



~2~ 



e 



(5.2) 

(5.3) 
(5.4) 
(5.5) 



(5.6) 



The eigenvectors of the preceding Hamiltonian operator corresponding to the eigen- 

2 



values A^ are the normalizable solutions of the differencial equation 



(5.7) 



We proved in Section 2 that if ip is a solution of (5/?), then the function z = will be 



a solution of the following Riccati equation 

dz 



(5.8) 



As it was stated in Section 3, such equation admits a nonlinear superposition principle 
based on the SL(2, M.) group, and therefore the general solution can be found by means of 
an appropriate factorization 



z(f) = exp(5f 2 L 2 ) exp(#iLi) exp(g L )(z) 



IC=o- 



(5.9) 



The functions go, gi, g2, are to be determined from the first order differential equation 

system 



go 

9i 

92 



e-x-gi 

-2g 



(5.10) 



together with the initial conditions #o(0) = gi(0) = #2(0) = 0. 
Let us first remak that the Riccati equation 

dz 



+ z — £ + \ = 0, 
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under the change of variables given by 

z = 2£v-Z i y = e, (5.H) 

becomes a new Riccati equation, 

dv n / 1 \ 1 — A , 

m^H*- 1 )-— - (5 ' 12) 

On the other side, the Riccati equation associated, according to the method described in 
Section 2, with the linear second order confluent hypergeometric differential equation 

yW" + (b-y) W'-aW = 0, (5.13) 

where a and b are constants and W and W" are the first and second derivative, respectively, 
of the function W(y), is 

dv 9 / b \ a , 

— + v 2 + v[--l)-- = 0. 5.14 

dy \y J y 

with 

V ~W' 



A simple comparison between ( {5.12j ) and (|5.14j) shows that both coincide when 



4 2 

It is well known that the general solution of ( |5.13| ) is given by 



1_A b=-. (5.15) 



W(V) = AM (l^, \, ij + By\M \, yj 

with A and B arbitrary constants, and M(a, b, y) is such that for large values of the variable 

y, 

M{a,b,y) = ie™ -a ( 6 + 1 ) (516) 

1 [b — a) 

+ ^\eyy a - b g(l-a,b-a,y), (5.17) 
r(a) 

where e = 1 if — n/2 < Argy < 37r/2 and e = — 1 when — 37r/2 < Argy < — 7r/2 and (7 
denotes the function 

= - a(a + l)-( fl + n-l)6(6 + l)-(6 + n-l) 

Therefore, the function behaves for large values of y as e y = unless that 

r(---| = ooandB = (5.19) 



4 4 / 

r(|-^J = 00 and A = 0, (5.20) 
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Now, we recall that Gamma function is such that T(0) = — oo and takes finite values 
when — z is not a positive integer number, while Y{—m) = oo for < m G Z. Therefore 
if we want to have a solution of (|5.13|) with values a and b satisfying (|5.15| ) being square 
integrable, it is necessary that 



1 _ A 

4 ~ 4 

3 _ A 

4 ~ 4 



= — m, and then A = 2n + 1 with n = 2m = 2, 4, 6, • • • (5.21) 
— m, and then A = 2n + 1 with n = 2m + 1 = 3, 5, 7, ■ ■ • , (5.22) 



or, 



- - - =0 and then A = 2n + 1 with n = (5.23) 
4 4 v ; 

3 A 

=0 and then A = 2n + 1 with n = 1. (5.24) 

4 4 v 7 

Consequently, a necessary condition for ( |5.13| ) with values a and 6 satisfying ( |5.15| ) to be 
square integrable is 

A = 2n+ 1, 

where n — 0, 1, 2, 3, 

Therefore, we will restrict ourselves to the case A = 2n+l. We will show that the solution 
for go can be obtained by induction on the index n. 

We first consider the case in which n is an even number. The first equation reduces for 
n = to 

go = e - 1 - gl (5.25) 

and then it is an easy matter to check that the solution we are looking for is go = — £■ 
In the case n = 2 the equation becomes 

9o = e - 5 - gl (5.26) 
and again it is easy to check that then the solution is 

9o=^^-Z- (5-27) 

In a similar way, when studying the cases n = 4, 6, ... we will see that for an even number 
n the solution satisfying go(0) = is given by 

*=f§H- (5 - 2s) 

with H n being the Hermite polynomial of order n and H' n means the derivative of H n with 
respect to £. In fact, using several properties of Hermite polynomials, we see that 

H'nH n — (H' n ) 2 — = (e-\)H 2 n -(H' n -H n 0\ (5.30) 
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and taking into account that H' n = 2nH n _i and the corresponding relation for the derivatives 
H" = 4n(n — l)H n ~2, we will get 



4n{n - l)H n „ 2 H n - H 2 n 



-\H% + 2£H n H' n , 
-XH n + 2£H' n . 



4n(n - l)H n - 2 - H n = 
Now, the recurrence relation H n = 2£H n _i — 2{n — l)_ff n _2, leads to 
2{n- l)# n _ 2 (2n + 1 - A) = 2£F n _ x (2n + 1 - A) , 



(5.31) 
(5.32) 



(5.33) 



and the right hand side vanishes because of A = 2n + 1. 

In summary, we have checked that for any even number n the solution such that go(0) = 
is given by 

= (5 ' 34) 



When introducing this value for g in the second equation ( 5.11 ) for gi we obtain the 
new equation 

g 1 = -2g = -2 f^-el, 



H„ 



which can be easily integrated 



9\ 



(5.35) 



(5.36) 



and as gi(0) = 0, C\ takes the value C\ = log(i/^(0)) = log £4. From the relation H n = 
2£H n ^i — 2(n — l)if n __2 we see that H n (0) = —2{n — l)i/ n _2(0) and iterating this reasoning 
we obtain the chain of relations 



# n _ 2 (0) = -2(ra - 3)tf n _ 4 (0) 
i2"„- 4 (0) = -2(n - 5)# n _ 6 (0) 



and therefore, 



and then, 



jf 2 (0) = -2(n - n + l)flo(O) = -2 • 1 

H n (0) = (-2(n-l)) (-2(n -3)) (-2(n-5))... (-2.1) 
= (-2) n/2 (n - 3) (n - 5) • • • 1, 

k 1 = H*{0) = [2 n ' 2 (n- l)(n-3)(n-5)---l 

2 n/2 n ( n ~ l)(7l-2)(7l-3)--- 

n(n-2)(n-4)---2 

2 



(5.37) 
(5.38) 
(5.39) 
(5.40) 



(5.41) 



nl 



n! 



(5.42) 
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So, the function gi is given by the expression 

9i = - log [hHle~^ 

Finally, the function g 2 is to be determined from the differential equation 

1 



92 



whose solution is 



92 



1 



where k 2 is to be chosen such that g 2 (0) = 0, i.e., 

1 



d£ + k 2 



k 7 



(5.43) 

(5.44) 
(5.45) 

(5.46) 



Putting together all previous results we get for the general solution z(£) the following 
expression: 



*(0 





=0) 


i - M 


(5=0) 




1 







+ 9o 

J7Z 



(C=o) 



1 + 



H -£ 



T2 



1 + 



We recall that the wave function ip was given by 
and from 



(5=o) 



(€=o) 



we obtain 



| «(eK = io g ir n -|-, 



(5.47) 



(5.48) 

(5.49) 
(5.50) 



Notice that when ip has a well defined parity, the quotient z — 4^ is an odd function 
and then the limit when £ — > of 2 should be zero. However, when ^ is a continuous odd 
function, then lim^o "0 = an d then the quotient cannot be finite. This was the reason 
for leaving aside the case A = In + 1 for an odd number n. In this last case it is convenient 



first to introduce the change of variable z = 1/v, and then the equation ( 5.8|) becomes 

| = 1 " {e - A) v\ (5.51) 
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and the corresponding system of differential equations 

g Q =l-(£»-A)sg 

91 =-2(e-*)go 

92 = -(£ 2 -A)e^. 

with initial conditions go(0) = gi(0) = (72(0) = 0. 
It is now easy to check that the solution for go is 



H„ 



9o 



that gi is given by 



and that g% is 



where 



g2 = - 



(e - A) h 



{H' n -iH n fe-e 
2n\ 



d£ + k 2 



(e - a) h 



Then, under the change z = 1/v the general solution 



(5.52) 



J «=o) 



1 - (02«) ( f=o) 

becomes after some computations the expresion given by ( |5.47| ) 



(5.53) 

(5.54) 
(5.55) 

(5.56) 
(5.57) 

(5.58) 



6 Conclusions 

The analysis developed in this paper of the reduction process leading to the Riccati equation 
starting from a second order linear differential equation just by a simple application of the 
general Lie theory of symmetry of differential equations, now a well established subject 
(see e.eg. Hl2|), allows us a better understanding of the so called factorization method, 
which is a method that have been used for finding new solvable potentials once one of such 
problem is known and that motivated the study of supersymmetric quantum mechanics. 
The Riccati equation was chosen here as the simplest example of first order differential 
equation systems admitting a nonlinear superposition principle and the deep reasons for the 
existence of such principle have been clarified in this paper. We hope that this will allow 
a geometric interpretation of the converse part of Lie-Scheffers theorem. This nonlinear 
superposition principle is very important and has played a very important role because 
it provides an explicit expression for the solution of nonlinear equations admitting such 
superposition principle. For a review see e.g. the review lecture given by Winternitz [jUjfl . 
We have also developed many examples of the use of Wei-Norman method of computation of 
the solution of differential equations in a Lie group and we have applied it as an academical 
example for reobtaining the eigenvalues and eigenvectors of the harmonic oscillator problem. 
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